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Abstract 

In supersymmetric standard models R-parity violating couplings are severely con- 
strained, since otherwise they would erase the existing baryon asymmetry before the 
electroweak transition. It is often claimed that this cosmological constraint can be 
circumvented if the baryon number and one of the lepton flavor numbers are sufficiently 
conserved in these R-parity violating couplings, because B/3 — Li for each lepton flavor 
is separately conserved by the sphaleron process. We discuss the effect of lepton flavor 
violation on the B — L conservation, and show that even tiny slepton mixing angles 
^12 ^ O{10^^) and 023, ^13 0(10^^) will spoil the separate B/3 — Li conservation. 
In particular, if lepton flavor violations are observed in experiments such as MEG and 
B-factories, it will imply that all the R-parity violating couplings must be suppressed to 
avoid the B — L erasure. We also discuss the implication for the decay of the lightest 
MSSM particle at the LHC. 



1 Introduction 



In most of supersymmetric (SUSY) standard models, R-parity is assumed to be exactly 

conserved, which prohibits the following baryon and lepton number violating operators, 

1 _ 1 

^RpV = -KjkLiLjEk + Xiji^LiQjDk + -X'ljfJJiDjDk + fiiLiHu ■ (1) 

R-parity guarantees the stability of the lightest SUSY particle (LSP) and plays a crucial role 
in LHC physics. In R-parity conserving models with the neutralino LSP, all the SUSY signals 
include a missing transverse momentum. If gravitino is the LSP, R-parity conservation leads 
to a long lifetime of the next-to-LSP (NLSP), which becomes a stable particle in collider 
scale. In fact, most of the LHC studies of SUSY signals crucially rely on these features. 

However, if R-parity is broken, the LSP is no longer stable, and the LHC signatures of 
SUSY events may drastically change. In particular, the decay length and the decay mode of 
the lightest ordinary SUSY particle (LOSP]0 crucially depend on the pattern and the size 
of the R-parity violating couplings [Ij. It is therefore very important to know what are the 
allowed R-parity violating couplings. 

The most stringent constraint on R-parity violating couplings comes from cosmology. 
Assuming that the baryon asymmetry is generated before the electroweak phase transition, 
baryon or lepton number violating processes induced by the R-parity breaking couplings 
(together with the sphaleron process [2]) would wash out the existing baryon asymmetry 
unless these couplings are sufficiently suppressed. The bound is roughly given by [3] 

A, A', A" < 0(10-^). (2) 

The bounds from the laboratory experiments and neutrino masses are much weaker than this 
cosmological bound (cf. [HH]). 

It is often argued, however, that the above cosmological constraint can be circumvented 
if the baryon number and one of the lepton flavor numbers are sufficiently conserved in these 
R-parity violating couplings, because B/3 — Li for each lepton flavor is separately conserved 
by the sphaleron process. (Here, B denotes the baryon number and is the lepton flavor 
number of the i-th generation.) For instance, if A'^^-^ < 10""^ is satisfied, then Ag^^, and X'^jf, 
can be much larger than 10~^. 

However, lepton flavor is not conserved in generic SUSY models because of the mixings in 

the slepton mass matrices. These slepton mixings can then erase the asymmetry between the 

^LOSP denotes the lightest SUSY particle among the superpartners of the standard model particles. If 
gravitino is the LSP, NLSP is the LOSP. If not, LSP is the LOSP. 
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lepton flavors, Lj — Lj, and hence {B — Li/ 3) — (B — Lj/3) [5j. In this paper, we investigate 
the effects of the lepton flavor violation (LFV) on the cosmological B — L conservation. We 
show that, if there exist slepton mixings as large as 9i2 ^ 10^^ and ^23(13) ^ 10^^, they would 
erase the lepton flavor asymmetry Lj — Lj and spoil the separate B — Li/ 3 conservation. It 
means that all the R-parity violating couplings must satisfy the cosmological constraint. In 
particular, this is the case if LFVs are observed at the current and future experiments such 
as MEG [3 and super B-factories [8j. We also reinvestigate the cosmological bounds on the 
R-parity violating couplings by solving the Boltzmann equations, and show that the trilinear 
couplings A, A', and A" must be smaller than 0(10^^ — 10~^), and the bilinear coupling should 
satisfy ^i/fi < 0(10^^). Finally, we discuss the implications of the bounds on the R-parity 
violation for the LHC phenomenology. 

2 Lepton flavor asymmetry in the early universe 

2.1 Slepton mixings and lepton flavor violation in the early universe 

In this section we discuss how fast the lepton flavor is mixed under LFV processes. Here, we 
do not introduce any R-parity violations to concentrate on the effect of the slepton mixings. 

Let us first consider the basis where the lepton Yukawa matrix is diagonal, which is 
familiar in the context of low energy LFV phenomenology. In this basis, we have a term 

W = hiLiEiHd (3) 

in the superpotential and soft slepton masses 

- £ = (m2 ),,L*Z, + (m|j,,f ;f,- (4) 

in the Lagrangian. To discuss the LFV effects in the early universe, however, it is more 
appropriate to take a basis where the slepton mass matrices are diagonal 0. Because the 
gaugino-slepton-lepton interactions are stronger than the Yukawa interactions, one should 
also diagonalize the gaugino interactions. Namely, leptons and sleptons are rotated by the 
same unitary matrices which diagonalize the slepton mass matrices. Assuming that the 
mixing angles are small, these rotations are expressed as 

Li = (UjJ.jL, U + Y^efjL,, (5) 
Ei = {U~,^),fij ^% + Y,e^!^j, (6) 
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where 6^-^ ~ ~^fi^ the mixing angles, and the hat denotes the new basis. Note that 
those mixing angles are different from the dimensionless parameters 

{5j^)ij = imfjij/{mf^)ii , {6ej,)ij = im^Jij/im1^)n , (7) 

which are familiar in the context of the LFV rare processes. In fact, the mixing angles O^- 
are enhanced compared to {6r (See discussion in Sec. 12. 2i ) 

In this new basis, the LFV effects appear only in the Yukawa couplings, which are given 



by 



Wlfy = Y,hijLiEjHd, (8) 



where hij = hiOf- + hjO^- . For instance, 

h23 = h2eg + hse^2 (9) 

(0.0061 + 0.10- 03^2) (^)- (10) 

We now estimate how much the lepton flavor asymmetry Lj — Lj is erased due to the 
above LFV interactions. To this end, we solve the Boltzmann equation for the evolution of 
Li — Lj. Here, for simplicity, we include only the effect of the higgsino decay and its inverse 
process, H «^ LiEj and H ^ LiEj, assuming that the higgsino is heavier than the sleptons. 
Other processes such as 2 — t- 2 scatterings and those with Higgs bosons may be comparably 
important, but it is expected that the bounds on the mixing angles will change only by 
order one factors. Note that these additional effects only strengthen the erasure effect, and 
therefore the bounds we will derive should be regarded as conservative ones. 

As is derived in Appendix IA.2[ the Boltzmann equation for the evolution of Li — Lj is 
given by [see Eq. (j60|) ] 

d _ 16(r,,+r,,) Fi(m^/T) ^ 
T^N,^-L, ^ F,{mj/T)+2^'^^-'^ ^^^^ 

where T is the temperature of the universe, H is the Hubble parameter, Fi{x) = x^Ki{x) 
with Ki{x) being the modified Bessel functions of the second kind. Nl--.^. is defined as 
Nii-Lj = {ni^ — nj^)/T^ — (n^^ — n-^.)/T^, where n^, and n-^_ denotes the lepton and 
anti-lepton number density in the i-th generation, respectively. The partial rate Fjj is given 

by 

r„. = Ltfm^fl-^^ , (12) 
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Figure 1: Time evolution of N^^-Ls for slepton mixing angles ^23'^ = 1 x 10~^, 3 x 10~^, and 
5 X 10~^, from the top to the bottom, for = 600 GeV, mj = 200 GeV, and tan/3 = 10. 
The vertical dashed line denotes the sphaleron decoupling temperature ~ 100 GeV. The 
normalization is arbitrary. The time evolution of N^-^-Ls for ^f/ ^ = (1 — 5) x 10~^ is almost 
the same. 

where and mj are the masses of higgsino and sleptons, respectively. We assume that the 
slepton masses are approximately the same. Note that the Boltzmann equation Eq. (jllh is 
symmetric under the exchange of the left-handed and right-handed slepton mixings, -f-)- 9^, 
i.e., they give the same effect on the evolution of Nl^-^.. 

In Fig. [H the time evolution of Nl^-l^ is shown for ^23'^ ~ (1 — 5) x 10~^, for nig = 600 
GeV, mj = 200 GeV, and tan (3 = 10. One can see that the flavor asymmetry is rapidly 
decreased for T < almost washed out for ^23'^ ^ 3 x 10^^. The time evolution of 

Ni^^L.^ for Oyl^ ~ (1 — 5) X 10^^ is essentially the same. 

In Fig. [2] and Fig. [3l we show the dilution factors 

as functions of the mixing angles 9- , where 100 GeV is the temperature when the 

sphaleron process is decoupled. In the numerical calculations, we take T^, = 100 GeV, = 
200, 600 and 1200 GeV, mj/rrif^ = 0.4 and 0.8, and tan/3 = 10. Note that the dilution 
effect is weaker for nig = 200 GeV than for 600 GeV. This is because for = 200 GeV the 
duration of the Lj — Lj erasure is shorter than for Trig = 600 GeV. 

One can see that the lepton flavor asymmetries L2 — -L3, Li — L3, and Li — L2 are washed 
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Lepton Mixing Angle 6 



Figure 2: The dilution factor D^j-Ls (-^ii-La) ^ ^ function of the slepton mixing angle 
(eifg) or (6'f3), for m^^ = 600, 200 and 1200 GeV, from the left to the right. The slepton 
mass mjis OAm^ for the solid lines and 0.8m ^ for the dashed lines. We took T^, = 100 GeV 
and tan /3 = 10. 

away for 

e^t > (0.3-1.0) X 10-^. (^)"\ (14) 

^ff > (0.3-1.0) X 10-^. (^)"\ (15) 

^ff > (0.6-1.6)xl0-^.(^)"\ (16) 

respectively. (We take the value where the dilution factor becomes D^.-l. ~ 0.01.) If any two 
of these inequalities are simultaneously satisfied, all lepton flavor numbers become essentially 
the same, Li = L2 = L3, and hence B — Li/3 = B — L2/3 = B — L3/3. As we will see in the 
next subsection, such slepton mixings are indeed naturally expected in generic SUSY models. 

2.2 Lepton Flavor Violation 

The lepton flavor violation is rather generic in a wide class of the SUSY models. In fact, the 
flavor structures of the neutrinos and quarks can induce mixtures of the slepton generations. 
In the see-saw models, the neutrino Yukawa interaction radiatively contributes to the left- 
handed slepton mass, while the right-handed one receives a correction from the CKM mixings 
above the GUT scale in the SUSY GUT models. By mediating the SUSY breaking effect at 
the Planck scale such as in the gravity mediation, the slepton mass matrices, thus, acquire 
the flavor mixing through the renormalization group evolution down to the weak scale, which 
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lO"*^ : , , i l < , , 

5x10-' 10-5 2x10-5 5x10-5 2x10-* 5x10-"* 

Lepton Mixing Angle 9 



Figure 3: The same as Fig. [2] but for the dilution factor Dl^-l^ as a function of the slepton 

r / o 

mixing angle 6-^2 ■ 

are approximately shown as 



jjt] - -^{Y*)ki{Yy)kj{2,ml + al)\n^^. 



("^D^. ^ -^Y^{VcKMMVcKM)U'^rnl + al)\n^, (17) 

where the ttt-q and ao are typical values of the scalar mass and the trilinear coupling. Here 
Mp, Mfi^ and Mh^ are the Planck scale, the mass of the k-ih. right-handed neutrino and the 
mass of the colored Higgs boson of the GUT, respectively. Obviously, the neutrino Yukawa 
couplings, Yy, and the CKM matrix, Vckm, induce the flavor violation even in the absence 
of the slepton mixings at the input scale. 

The flavor mixings in Eq. pzp are represented in terms of the parameters at the weak 
scale and those of the right-handed neutrino. It is known that the neutrino mass is expressed 



as iYu)ij = 1 / {Hu) \J M/j . y/'^i'k ^jk ■^ith the vacuum expectation value of the up-type 
Higgs boson, (Hu), the light neutrino mass, mi,^, and the PMNS matrix, Ujk [6]. Here, Ujk 
and Rij satisfying Rn^-Rj/^ = 6ij generally have a nontrivial flavor structure. On the other 
hand, we assumed the minimal Yukawa interactions above the GUT scale in the latter line of 
Eq. pT|) . Besides, we took only the top Yukawa coupling, Yt, into the account because the 
up-type quark masses are hierarchical. Noting that the diagonal components of the slepton 
mass matrices are approximated to be ijiq, we estimate the flavor mixings as 

Mr 



-5 



IQio GeV 



{Sen)2i ~ 3 X 10-^ {6^^)32 ^ 3 X 10-2, {5^^)31 ~ 7 X 10-3, (18) 
where {SA)ij is defined as {^Ajij = (™'A)«i/("^yi)«« 7^ i)' have assumed Oq = uiq and 
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Mh, ^2x 10^^ GeV. We also chose the typical mass of the light neutrino to be O.leV and 
Rij = 0(1) discarding cancellations. 

We find that the typical size of the flavor violation is large enough to erase the flavor 
dependence of B/3 — Li. Namely, we cannot evade the wash-out condition by postulating 
the flavorful R-parity violation. This conclusion may be seen simply by equating 9-j in 
Eqs. (fn|) " p^ with {^Jj^i^^)ij in Eq. ([TH]) . Actually, the former becomes much larger than 
the latter when the slepton (sneutrino) is almost degenerate. They have a relation as 

where Am\ is the difference of the slepton (sneutrino) mass eigenvalues, while rn- and 
in the numerators denote the average of the mass eigenvalues. In many models, the slepton 
(sneutrino) mass matrices are set to be universal at high scale to avoid too large lepton- 
flavor violations. Then, the mass difference stems from the Yukawa interactions through the 
renormalization group evolution. Since the effect is very small, 6^-^ becomes much larger 
than ((5a) jj) leading to the rapid mixture of the sleptons in the early universe. 

The above flavor mixings in Eq. (jl7p induce lepton-flavor violating decays of muon and 
tau-lepton. The branching ratio of the /i — )• 67 decay is approximately obtained as |9l[^ 

Br(^ ^ e7) ~ 10-(--3) I ^^}A V tan^V f^I^\ ^ (20) 
" y 10-4 j \ m ) V400GeV/ ^ ' 

where mgoft is a typical mass of the sleptons, sneutrinos, neutralinos and charginos. Here, 
the coefficient is quite sensitive to the details of the sparticle mass spectrum. In the presence 
of the right-handed slepton mixing, ((^6^)21, we checked that the additional contribution is 
smaller by a factor or by an order of magnitude than Eq. (|2U|) . We notice that the decay 
rate is proportional to ((5^)^^ rather than {Oij)"^. Thus, when Am^ is small, the rate becomes 
suppressed for fixed 6ij, as is expected from the GIM mechanism. 

In Sec. 12.11 we found that the flavor mixing between Li and L2 becomes effective in 
the early universe when satisfles Eq. (jl6p . From Eq. (jl9p and Eq. ()20p . this condition 

corresponds to 

or similar condition with fn-^. On the other hand, the current experimental bound and the 
^We neglect the contributions from multiple flavor mixings such as ((SA)23(5yi)3i. 
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future sensitivity are 

Br(^^e7) < 1.2 x 10"i\ 

Brifi^ej) > 10-'^ (22) 

from the MEGA [TO] and MEG [7] experiments, respectively. Comparing these two results, 
we conclude that if we measure — )• 67 decay in future, it suggests that the mixture of the 
lepton flavor takes place effectively in the early universe. 

Let us next consider the flavor violating decays of the tau lepton such as r — t- /u(e)7. The 
branching ratios are estimated as |9j 

B,. . M.).) ^ 10-..3-.. {^^)\^y i^y. (23) 

The right-handed slepton mixing contributes again by a factor or an order of magnitude 
smaller than the left-handed one. On the other hand, the current experimental bounds are 

Br(r fi'j) < 4.5 x 10"^ 

Br(T ^ 67) < 1.1 X 10"^ (24) 

from Belle and BaBar [11]. The sensitivity of the branching ratio is expected to be improved 
by an order of magnitude in future super B-factories |8j. Compared with the results of 
Eqs. ()14p and (jlSp . we find that the cosmological flavor mixing of the tau lepton becomes 
effective when the branching ratios satisfy 

B,.(.^M.)7)>10-<-«)(^)^(^)-'. (25) 

This lower bound is much smaller than the future sensitivity. Therefore, if we observe the 
tau- lepton flavor violation in future colliders, then the conditions Eqs. p4p and/or (llSp 
are satisfied, and the tau-lepton fiavor number is not independently conserved in the early 
universe. 

3 Implications for the R-parity violation 

As we have shown in Sec. 12. H if the slepton mixing angles satisfy at least two of Eqs. (I14p - 
(I16p . then all the lepton flavor asymmetries are equilibrated, i.e., Li = L2 = L3. In fact, 
as was discussed in Sec. 12.21 sizable slepton mixings are expected in a wide class of SUSY 
models, which are sufficiently large to satisfy Eqs. (jl4p - (|16p . In this section, we discuss the 
bounds on the R-parity violating couplings in the presence of such lepton flavor violations. 
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3.1 Cosmological bounds on the R-parity violation in the presence of 
slepton mixings 



We assume that at least two of Eqs. (ll4p -( fT6]) are satisfied, and hence all B — Li/3 are 
equilibrated. Then, in order to avoid the baryon erasure, any of B — Li/3 violating processes 
should not become effective before the electroweak transition. 



We calculate the dilution factor 

Db-l 



(26) 



iVB_L(r»r,) 

as functions of the R-parity violating couplings Ajjfc, X'-f,, A-'^-^, and /ijj^ The corresponding 
Boltzmann equations are shown in Appendix lA. 31 The results are shown in Figs. SHZl Here, 
for simplicity, we have assumed that all sleptons and all squarks have the same masses mj 
and TTiq, respectively. 

From the figures, one can see that the couplings should satisfy 

< (4-5)xl0-', 



ijk\ 



ijk 



IJ2Kk\ 

ijk 
ijk 



E 



< 



< 



< 



-7 



(3 - 6) X 10 



(0.6 - 1) X 10" 



(1 -2) X 10" 



tan /3 
10 



(27) 
(28) 
(29) 

(30) 



for niq ~ 



200 - 1200 GeV and nij ~ 100 - 400 GeV. (Again, we took the value where the 
dilution of the B — L becomes Db-l — 0.01.) We should note that the bound on the U D D 
coupling X"jf^ in Eq. (i27|) applies even without the lepton flavor violation. 



3.2 Imphcations for coHider phenomenology 



In the presence of slepton mixings, all the R-parity violating couplings must satisfy Eqs. (|2< 
(j30p in order to avoid the baryon erasure. Interestingly, this means that the LOSP has a long 
decay length at the LHC. For instance, suppose that the LOSP is the stau, mainly consisting 
of the right-handed stau. If the LLE coupling Ajjs saturates the cosmological bound Eq. ([29j) . 



^We do not discuss the bounds on the R-parity violating soft terms, for simplicity. 
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5x10"* 10-' 2x10-' 5x10-' 10-'' 

U DD coupling A" 



Figure 4: The dilution factor Db-l in the presence of an R-parity violating term \"U iU jD^ 
for rriq = 600, 200 and 1200 GeV, from the left to the right. We took mj = lOOGeV and 
= SOOGeV, but this result is nearly independent of these masses. Other parameters are: 
tan /3 = 10 and = lOOGeV. 

the decay length of the stau becomes 

„,.50.,n(^)-^(^)-. (31) 

This is comparable to the tau-lepton decay length (cr,- ~ 87 /^m), which can be probed at 
the LHC. Note that this is the shortest possible decay length, and in general a (much) longer 
decay length is expected. If the dominant decay of the stau is caused by Ajjfc {k ^ 3) or the 
LQD coupling X'-jf^, the decay length becomes longer since the decay rate is suppressed by 
the left-right mixing of the stau and/or the flavor mixing. 

Similar results are obtained for other LOSP cases. Depending on what is the LOSP and 
the pattern of the R-parity breaking, the dominant decay mode of the LOSP can become 
three- or four-body decay [Ij, which makes the decay length even longer. It is important to 
study the LHC phenomenology of R-parity violating SUSY models under the cosmological 
bounds Eqs. (p7|) - ([30]) for different LOSP candidates and the different pattern of R-parity 
violating couplings, and we leave it for future work. 

4 Conclusion and discussions 

In this paper, implications of the slepton flavor mixings for the cosmological constraints on 
the R-parity violating couplings are discussed. We have shown that tiny slepton mixing angles 
6*12 ^ 0(10"^) and 023,^13 ^ 0(10^^) will spoil the separate B/3 — Li conservation. Such 
slepton mixings are indeed expected in generic SUSY models. On the other hand, if lepton 
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5x10"* 10-' 2xUr' 5x10-' 10-'' 

LQD coupling A' 

Figure 5: The same as Fig. [J] for X' LiQjD]^ interaction. Parameters are the same as Fig. HI 
mj and mjr hardly affect the result again. 




LLE coupling A 

Figure 6: The same as Fig. H] for XLiLjEk interaction, mj = 400, 200 and 100 GeV, from 
the left to the right, rrig = 600GeV and = SOOGeV. In this case the result depends on 
mj, and is almost independent of the other masses. 

flavor violations are observed in experiments such as MEG and B-factories, it will imply that 
all the R-parity violating couplings must be suppressed to avoid the B — L erasure. 

We also reinvestigated the cosmological constraints on the R-parity violating couplings 
in the presence of lepton flavor mixings in the slepton sector, and showed that the R-parity 
violating couplings must satisfy Eqs. (p7|) - ([30]) . It then suggests that the LSP (or NLSP in 
the gravitino LSP scenario) has a long decay length, which can be probed at the LHC. 

Interestingly, for such a small R-parity violation, the gravitino becomes an attractive 
candidate for the dark matter [12^ I13j. where the gravitino lifetime can be much longer than 
the age of the universe due to the double suppression by the inverse Planck mass and the 



12 



10"' 2x10"' 



5x10"' 10"'' IxlO"*" 
LH coupling e 



Figure 7: The same as Fig. S] in the presence of a bilinear R-parity violating term //jLjff^ as 
a function of = ^i/ ^. The masses and the other parameters are the same as Fig. [H The 
result is nearly independent of mj, rng, and the generation index i. 

R-parity breaking coupling. 

So far in this paper, we have assumed that the baryon asymmetry was generated before 
the electroweak transition, and found that if LFV events are observed, then the R-parity 
violation must be so small that the LSP becomes long-lived in the LHC. This argument can 
be put the other way around. If LFVs are discovered and also a sizable R-parity violation 
is observed at the LHC, then it means that the baryon asymmetry of the universe must 
have been generated after the electroweak transition, as in the electroweak baryogenesis and 
in Affleck-Dine baryogenesis with a long lived condensate or Q-ball. This is also quite an 
exciting possibility. 
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A Boltzmann equations for lepton flavor and R-parity violat- 
ing processes 



In this Appendix we derive the Boltzmann equations for lepton flavor and R-parity violating 
processes in the MSSM. 

A.l Boltzmann equation 

The time evolution of the number density nA_ of certain particle A obeys the Boltzmann 
equation. When we consider only the part of the time evolution induced by a process X ^ 
AY, that is, the decay of some particle X and its inverse process, the Boltzmann equation is 



at 



nx (Tx^ay ^ - riAUY ^ (cxu) ^y^x ^ , (32) 



X^AY 

where H is the Hubble parameter, and ( ) denotes thermal average. We assume that X, A, 
and Y are all in thermal bath, and discuss the effect of (very weak) X ^ AY process. Using 
the "yield" N = n/T^ as a variable, Eq. ([32]) becomeJll 

1 



X^AY '- 



nxiTx^AY) - nAnY{{crv)AY^x) ■ (33) 



For simplicity, we approximate the distributions of the particles as the Maxwell-Boltzmann 
distributions Ifl Then, the number density n of a particle is 



where //, g and m are the chemical potential, the degree of freedom, and the mass of the 
particle, respectively, and Fi{x) is defined as 

Fi{x) = x'K.ix) (35) 

through the modified Bessel function Ki{x) of the second kind. For a massless particle, 
-^2(0) = 2. The rate of the decay process is 

'^^Px f Ex- H.x\ mx ^0 



nx{Tx^AY / =9x J exp J -^'Tx^ay (36) 

-s^rU-«~F,{!!^)e.,{l^). (37) 



*We have used dT/dt = ~HT, assuming for simplicity that the effective degrees of freedom is 
constant. 

^Even if one uses the Bose-Einstein and the Fermi-Dirac distributions, the foUowing discussion is almost 
unchanged. In this case, one has to include the statistical factors for the final states, and the functions Fi(x) 
are replaced with modified functions. 
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Here Tj^^^^y is the partial decay rate in X's rest frame. Similarly one can calculate the rate 
of the inverse decay process, which is 

T3 



(38) 



Thus when the process X AY is in thermal equilibrium, which means nx = fJ^A + fJ-Y, the 
decay process and inverse decay process occurs with the same rate. 

Here we introduce the difference between decay and inverse decay process, or "effective 
decay rate," as follows: 



X ^ AY)) = nxiTx^AY) - "A^y ( (o-i')Ay^x 



9xr'x^AY^Fi{^) 



The equation Eq. (p3|) is now 



1 



exp j - exp I 



X ^AY 



(39) 
(40) 

(41) 

X^AY 

Moreover, as the chemical potential /x^ of the antiparticle is equal to —^a, the effective rate 
of the processes of antiparticles are 

r3 



exp 



-^J'X 



exp 



fJ-A + fJ-Y 



and therefore 



X^AY 



gxT^ 



Xt± AY 



X^AY 



Xt± AY 



sinh (^^^^ ~ sinh 



MA + [i-Y 



(42) 

(43) 

(44) 



A. 2 Lepton flavor violation 

Here, we derive the Boltzmann equation for the LFV process in the early universe. As 
an example, we consider as LFV processes those induced by the following term in the 
superpotential: 

W = h2sL2E3Hd . (45) 

For simplicity, we discuss only the decays and inverse decays of the higgsinos H, and assume 
that all sleptons have the same mass mj (< m^). 
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We define the asymmetry N^j^-^ of a supermultiplet as, e.g. for ^ul, 

= {N,, - Nj^) + (iV^ - AT^. ) . (46) 
Since leptons are massless before the electroweak transition, the asymmetry is 



iVi 



[ml] 27r2 
1 

~2 



exp 



T 
sinh 



exp 



+ 



ML 



27r2 



exp 



\ T 



exp 



V T 



(47) 



and its time evolution induced by LFV processes is described as 



1 



LFV 



-(their antiparticles' processes) 



where 



2 r 
-2-— T7^i 



r 



sinh ( - 



-ms 1 



m-f 



sinh 



+ ^^E■, 



(48) 
(49) 

(50) 



327r 

is the partial decay rate of each process, which is the same for all four processes. Similarly 



LFV 



LFV 



2 r 

'2 ■ ~2 -^1 \ ^ 



, 2 r ^ fnijr 
-4 • F^' ^ 



sinh 
sinh 



T 



— sinh 
sinh 



^^L2 + 
T 

T 



(51) 
(52) 



Now let us consider the asymmetry of each lepton flavor, which is defined as 
and so on. Prom Eq. ()47p . they are given by 



Ni 



+ 2 



(53) 



(54) 



Under the LFV interaction L2E3H11, the time evolution of the difference N2 — is given by 



T—{N2 - Ns) = — + iV[.^] - iV[^]) 



16 r 



sinh 



+ sinh 



16 r 



f^Hg + ^^2 + /^g3 

T 



(55) 
(56) 

(57) 
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where we have used <C T. On the other hand, reactions mediated by the diagonal lepton 
Yukawa couphngs are in thermal equilibrium for T < 10^ GeV, which leads to 



and hence 



16 r 



TT^H 



3r 



Therefore, from Eq. (|54p one obtains 



(58) 



(59) 



(60) 



A. 3 R-parity violation 

The evolution of the B — L asymmetry under the R-parity violating interactions can be 
discussed in the similar way as Sec. IA.2I To this end, we first discuss the relations between 
the chemical potentials and the B — L asymmetry in Sec. IA.3.11 We then discuss the time 
evolution of the B—L asymmetry in the presence of R-parity violating couplings in Sec. lA.3."2l - 
Sec. ESl 



A. 3.1 B — L asymmetry and chemical potentials 

For 10^ GeV < T < 10^ GeV, the reactions mediated by Yukawa and gauge interactions as 
well as the sphaleron process are all in thermal equilibrium, leading to the following relations 
between the chemical potentials of the MSSM particles 



+ ^IJ'L, = . 



(61) 
(62) 
(63) 
(64) 



Note that the chemical potentials of gauge bosons and gauginos vanish, and hence particles 
in the same super- and gauge multiplet have the same chemical potentials. Moreover, quark 
mixings make their chemical potentials independent of the flavors. Thus, all the chemical 
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potentials can be expressed in terms of four of them, e.g., //l- and hh^'- 

i 

i 



I 



f^E, = -f^L, - f^Ha ■ (68) 

Let us define the asymmetry of a supermultiplet N^y^^ as 

N^^^^{Na-N^) + {N^-N^.). (69) 



Since quarks and leptons are massless before the electroweak transition, from Eq. p4p one 
obtains 



Ni 



• ffeff (^) sinh (^) , 5cff (x) = 2 + F2(x) , (70) 

for A = Q,U,D,Li,Ei. Their degrees of freedom are given by (gg, gij, 9^, gL,, gE~) = 
(18,9,9,2,1). For simplicity, we assume that all the squarks and the sleptons have the 
same masses rrig and mj, respectively. Then, the hypercharge conservation 

Y^YaN^a]=0 (71) 



leads to 



9eS[-jrj I J, ) + ffeff y-jT 



'mrr 

5effl 



(^)p^ =0, (72) 



where we have used fi <^ T and neglected the Higgs boson masses for simplicity. Therefore, 
from Eqs. ([65]) - ([68]) we obtain 



t^Ha = -Ch^ iT) ^ /XL, , (73) 

i 

where 

.^,1 2geff(mg-/T) + 6geff(mj/T) 

""'^ ' 3' 9g,s{m^/T) + 3g,s{mj/T) + 2g,s{m^/T)- ^ > 

Together with Eqs. ()65p - (|68p . all the chemical potentials are now given in terms of /i/,. . Note 
that these relations are independent of the existence of lepton flavor and R-parity violations. 
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We define the B — L asymmetry as Nb-l = {ns — riLj/T^ ■ From Eq. ([70]) . it is given by, 
for ij, <^T, 



where 



-^Cb-l(T)^^, (75) 



4 / 777/ — \ / \ / Tft — \ 

Cb-l{T) = 3 5eff(^) + 3(1 - CH.iT)) geff(y) (76) 
and CHaiT) is given by Eq. (fTI]) . 



A. 3. 2 U DD interaction 

In the presence of the superpotential 



W = h'l^,U,D,Dk, (77) 



the time evolution of Nb-l is given by 



where 



T^{Nb-l) ^ Fi^— j (78) 

- ^ Vr F 3 + 9Cj,,(r) 

%:D.D.=I^I^^'.'^lK- (80) 

We should emphasize that Eq. ()79p holds even in the absence of the lepton flavor violation. 
A. 3. 3 LLE interaction 

Here we assume that the lepton flavor asymmetries vanish because of the lepton flavor 
violation, that is, we use = = m addition to Eq. (f75]l . Under this assumption, 
the time evolution of Nb-l under the superpotential 

W = ^XijkLiLjEk (81) 

is described as 

T±Nb-l ^ _i_i5!£^^ilMii.,(!^~)^/f£±ifl (82) 

dT 7r2 H ) T ^ ' 

^ F C^A i + 3g//.(r) ^ 



ijk 
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where 

r, 

A. 3. 4 LQD interaction 



^L,L,E, = ^\^ijk\^mi- (84) 



Here also we assume the vanishment of lepton flavor asymmetries. The time evolution of 
Nb-l under the superpotential 

W = K^kL^QiDk (85) 



IS 



'l^^^-^ 2 

dT -K^ 



H J H \T 



H 

ijk 



where 



A. 3. 5 Bilinear R-parity violation 

The bilinear R-parity violating term fiiLiHu induces, through the Li-H^ mixings, effective 
trilinear couplings Xijk and A^^-^. Then, the time evolution oi B — L can be discussed by using 
the Boltzmann equations in Sec. IA.3.31 and Sec. IA.3.41 
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